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Abstract 

The gauge dependence of the renormalization group functions of the 
Ginzburg-Landau model is investigated. The analysis is done by means of 
the Ward-Takahashi identities. After defining the superconducting order pa- 
rameter, it is shown that its exponent (3 is in fact gauge independent. This 
happens because in d = 3 the Landau gauge is the only gauge having a phys- 
ical meaning, a property not shared by the four-dimensional model where 
any gauge choice is possible. The analysis is done in both the context of the 
e-expansion and in the fixed dimension approach. It is pointed out the differ- 
ences that arise in both of these approaches concerning the gauge dependence. 
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I. INTRODUCTION 



It is often stated in the literature that the superconducting order parameter in a 
Ginzburg-Landau (GL) model has no physical meaning because it is a non-gauge invari- 
ant quantity. This statement is obvious if by order parameter we mean < 4>(x) >, the 
expectation value of the complex scalar field in the GL model. Since physical quantities 
should be gauge invariant, < <f)(x) > cannot be considered as a physical quantity. It is 
desirable, however, to have a gauge invariant definition of order parameter which character- 
izes the superconducting phase transition. One such definition exists already and has been 
proposed by Kennedy and King ]I[ in the context of the lattice superconductor. The order 
parameter proposed by them is given by G = limix.yi^oo < (f)^(f)* >, where X is a lattice 
field. This order parameter describes a phase transition in the lattice GL model. However, 
it is not clear that this phase transition coincides exactly with the normal-superconducting 
transition or, in the language of particle physics, the Higgs transition. 

In this paper we will discuss the physical meaning of the superconducting order parameter 
in the context of a continuum GL model. The aim of this paper is to discuss in a deeper 
way the questions addressed by one of us in a recent letter [Q]. It will be shown that a 
gauge invariant definition of order parameter, consistent with the traditional definition of 
the critical exponent /3, is in fact possible. Less obvious is the gauge independence of (3 in 
different schemes of renormalization. For instance, the gauge dependence obtained through 
dimensional continuation, the e-expansion, is very different from the one obtained using 
a fixed dimension approach. This approach, though less controlled, gives better values to 
the critical exponents. In fact, the renomalization group (RG) calculations based on the e- 
expansion does not give very good results for the superconducting phase transition, specially 
in what concerns the type II regime. Indeed, this regime is not readily accessible by the 
e-expansion which predicts a weak first order transition whatever the value of the Ginzburg 
parameter k 0. As shown by Dasgupta and Halperin Q using duality arguments in a lattice 
model, the weak first order scenario does not hold in the type II regime where the transition 
should be expected to be second order. In order to describe correctly the critical behavior of 
superconductors using the e-expansion, it is necessary to use resummations techniques such 
as the Pade-Borel resummation used by Folk and Holovatch ||. 

The plan of this paper is as follows. In section II we will discuss the RG equations 
for the GL model and define the order parameter. In section III we will use the Ward- 
Takahashi (WT) identities to establish the gauge dependence of the order parameter exactly. 
It will be shown that in a fixed dimension d = 3 RG only the Landau gauge a = has a 
physical meaning. This result will be anticipated already in section II by looking the 1-loop 
approximation to the RG functions. However, it will only be proved in section III to all 
orders in perturbation theory. Section IV concludes the paper. 



II. RENORMALIZATION GROUP IN THE GL MODEL 



Let us consider the following bare action for the GL model: 




d d x \\f{ + P»t(D> ) + ^LAlAl + ml\M 2 + f I0o| 4 | + S gft (11.1) 
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where the zeroes denote bare quantities, F 2 is a short for Fq U Fq U and D° = + ieoA ^. 
The S g f is the gauge fixing part and is given by 



(II.2) 



We have added a mass to the gauge field in order to regularize the infrared divergences 
arising from the gauge field propagator. This mass breaks explicitly the gauge invariance 
(also broken by S g f), which is recovered at the critical point. Adding a mass to the gauge 
field does not affect the renormalizability of the model. This should be contrasted with the 
non-abelian case, where the only way to provide a mass to the gauge field without destroying 
renormalizability is through the Higgs mechanism ||. 

The renormalized action is obtained from the bare action by rewriting it in terms of 
renormalized quantities. This is achieved through the introduction of renormalization con- 
stants. We introduce the field renormalization in a standard way ||, — Z^' 2 A° and 

<p = 1//2 0o- The renormalized action including the other renormalization constants is 
given by 



S = J d d x 
+ Z m m 



^F 2 + Z^D^)\D, 



Z M M 2 



-A 2 



Z„ u , 



2| ^ 2 + ^l0| 4 + §(^) 2 



(II.3) 



The renormalization constants Zm and Z a above are in fact superflous. Indeed, we have 
that Zm = Z a = 1 which means that the corresponding terms in the action, the gauge field 
mass term and the gauge fixing term, are not renormalized. This result is easily checked 
by studying the WT identities as follows. By adding sources to the corresponding fields we 
obtain the following WT identity: 



M 2 A ) d, 
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SJ^x) 



+ ie 



x 



6Ji(x) 



J(x) 
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SJ(x) 



\w(J tt ,J,J) = dMx), (II.4) 



where W = log Z, Z being the generating functional of correlation functions defined by 



Z{Jp, J\ J) = J D^DtpDAfj, exp -S + J d d x{J^ + jty + $J) 



(11.5) 



When the sources are zero Z gives the partition function. W is the generating functional of 
the connected correlation functions. The Legendre transform of V of W is given by 



W{J^ J\ J) + rV, <p, a M ) = J d d x(J^ + jV + <P*J) 



(II.6) 



where the functional effective action r(</p', tp, a M ) is the generator of the 1-particle irreducible 
functions. It satisfies a WT identity which is the Legendre transform of ( |11.4j ): 



A-M 2 )d fl a fl (x) + d, 



5T 



+ ie 



o. 



(II.7) 
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Both WT identities ( [II. 4j ) and ( |II.7| ) are linear. Thus, if we perform a loop expansion, 



T = Ysi^h where Ti is the /-loop correction to the effective action, it follows that each Ti 
satisfy separately a WT identity ( |L1.7| ). The zero- loop contribution T is the classical action 
S, which satisfies the inhomogeneous WT identity (|II.7|) , that is, the identity containing the 
term (a _1 A — M 2 )d fl a /I (x). The loop corrections, on the other hand, satisfy the homogeneous 
WT identity. The non-gauge invariant terms are all included in the inhomogeous WT 
identity, do not receive any correction due to the fluctuations and are in this way not 
renormalized, that is, the coresponding counterterms are zero. The only non-gauge invariant 
terms in the action are the gauge field mass term and the gauge fixing term. As a consequence 
of these reasonings we obtain that Zm — Z a — 1. This implies that M 2 



Z A x a§. From ( [II. 3|) we deduce also m 2 = Z^Z^m^, e 2 = Z A el and u = Z^Z~ U 



Z A M 2 and 
l u . 



It is useful at this point to fix the renormalization conditions necessary to define the 
corresponding renormalized quantities: 



li 2) (0) 



m 



i j ■ 



(II.8) 



dT 



(2) 



dp 2 



'131 



p 2 =0 



{11.9) 



r i 4 L(° ; °> °> °) = ( s ij$ki + dikfiji + Sii6 jk )u, 



T^(O) = M 2 5 



(11.10) 

(11.11) 



_d_ 

dp 2 



p2 



F(2) 
jxv 



d-h 



(11.12) 



d p^ 
dp 2 p 2 



F(2) 



p 2 =0 



(11.13) 



where the latin indices represent the components of the scalar field. 

We are interested in the behavior when m — > 0, the infrared behavior. Thus, we assume 
that ml oc t, t = (T — T c )/T c being the reduced temperature, and m = £ _1 ~ t v as t — > 0, 
£ being the correlation lenght. The RG equations will be obtained by differentiating with 
respect to lnm with all bare quantities fixed, except for m\. Let us define the dimensionless 



couplings / = m° 
readily obtained: 



-v, 



u 



m d A u and v = m/M. The following exact flow equations are 



m- 



dM 2 

dm 
da 

dm 
df_ 

dm 
dv 

dm 



m 



m 



m 



VaM 2 , 



-r] A a, 



(VA + d-Qf, 



1 T lA\ 

T 



(11.14) 
(11.15) 
(11.16) 
(11.17) 
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where we have introduced the RG function t]a which at the fixed point gives the anomalous 
dimension of the gauge field. It is defined by 

VA = m^-\ogZ A . (11.18) 
At 1-loop order and in fixed dimension d = 3we obtain 

i* = h < IU9) 

From Eq. (|II.16|) we obtain that a charged fixed point should correspond to t]a — 4 — d. 
In d — 4 this corresponds necessarily to zero charge while in d = 3 we have /# = 247r, a 
non-zero charge. The only fixed point to Eq. ( [[I.15| ) corresponding to a charged fixed point 
in d = 3 is a* = 0, that is, the Landau gauge. Due to the negative sign in (|II.15|) , any a / 



in the flow diagram will flow to infinity. Fig. 1 shows the flow diagram in the /a-plane in 
the 1-loop approximation in d = 3. Note that the line a = contains an attractive flow 
towards the fixed point = 247r. Any a ^ 0, no matter how small, will flow away. This 
does not happen at d = 4 where f* = and a = a* arbitrary is a line of fixed points. This 
1-loop calculation suggests that in d = 3 the Landau gauge is the only physical gauge. 
The 1-loop beta function for the coupling u is given by 

(fit i) 

(3u = m— = (2 V<t> - l)u + ^-u 2 + —f 2 . (11.20) 

Om 07T 27T 

This beta function does not exhibit any gauge dependence, which is cancelled out between 
the rj^u and f 2 term. In ( |11.20| ) is defined as m<91n Z^/dm, which computed in the Landau 
gauge gives 

** = -S(ITJF (IL21) 

We will see in section III that the gauge independence of /3# is valid to all orders, though 
does depend on the gauge. 

The flow equations we have obtained are similar to the ones obtained by Herbut in 
the context of the continuum dual GL model PHT2"|j. The difference is that in our case the 
charge e is not meant to be the dual charge and, as a consequence, our mass M has nothing 
to do with the photon mass (anyway, this is a subtle point, even in the dual approach; see 
ref. 0). 

Let us discuss the critical behavior that arises from the so defined GL model. From an 
experimental point of view, presently only the exponents a (specific heat) and v (correlation 
lenght) are accessible. The order parameter exponent f3 seems not to be directly accessible. 
This is the main point we would like to discuss. At this point, an important remark is 
in order. The critical exponents are obtained from the singular behavior of the effective 
action (the free energy) and correlation functions. The WT identities imply that the singular 
behavior is the same, whatever t < (the broken symmetry regime or ordered phase) or t > 
(symmetric regime or disordered phase). This means that the exponents have the same value 
below or above T c . This remark is very important concerning the superconductors. The point 
is that below the transition the photon becomes massive, phenomenon known in particle 
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physics as the Higgs mechanism 0. This mass is generated spontaneously and not added 
by hand as our mass M. Since the singular behavior is the same as in the symmetric phase, 
the renormalization proceeds exactly as in the symmetric phase. Thus, this mechanism is 
very important in the context of non-abelian gauge theories, where adding a mass by hand 
destroys renormalizability. The photon mass generated by the Higgs mechanism in the GL 
model corresponds to the inverse of the penetration depth A. This lenght is known to scales 
as the correlation lenght £ J7J, IB| , Yl ] , that is, A ~ £ ~ \t\ u . A is a quantity that arises only for 
t < but its exponent, being the same appearing in the scaling of £, can be evaluated for 
t > or even for t = (the critical point), this last case needing renormalization conditions 
different from the ones we use here due to the infrared divergences [T3|]T~| • 

The v exponent is obtained easily by considering the flow of Wq (remember that tuq is not 



kept fixed under RG). From the definition of its renormalized counterpart m 2 
we obtain 



m *f = (2 + "" 



where 



dlnZ n 



m- 



dm 



[11.22) 



(11.23) 



Eq. fL22D implies that near the phase transition 



m 



0) wig ~ m 2+r,m v , where r)^ and 77 



are the fixed point values of r/ ri 
we obtain immediately 



and 7?0, respectively. Since by definition cx t and m ~ t" 



(11.24) 



The v exponent can be accurately measured through a direct measurement of the pene- 
tration depth |15|]. In the accessible critical region, the fluctuations are governed by the 3D 
XY fixed point |16|], which is a neutral fixed point. Then, it results that p s ~ A -2 ~ t 2u , 
where p s is the superfluid density |12|,|16j and we have defined the exponent v' of the pen- 
etration depth. Using the Josephson relation p s ~ [ 1_S | for the d — 3, we obtain 
v' = vj2. The exponent v' obtained experimentally in bulk samples of Y Ba2Cu^Ois is 



v' ~ 1/3 fl5| and therefore we obtain v pa 2/3, corresponding in this way to a 3D XY uni- 
versality class. This very same value of v is expected for the charged transition, the so called 
"inverted" 3D XY universality class |4j]. In fact, this is confirmed by recent RG calculations 

171 



3j] and Montecarlo simulations 
present framework. We have at 1-loop order that 



Let us evaluate v from Eq. ( |11.24[ ) in the 



Vn 



47T* 



(11.25) 



The infrared stable fixed point corresponds to /* = 247r, = 871/ '5 and v* = 0. This gives 
1] = and v = 0.63 at 1-loop order. This result can be systematically improved. Indeed, the 
fixed point will always have a 3D XY value because every power of / will be multiplied by 
a function of v which goes to zero as m — > 0, since v ~ m 1//2 . Therefore, we obtain already at 
2-loops v pa 2/3. It should be noted that this behavior is obtained because from Eq. ([II. 14|) 
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we have that M 2 ~ m 4 ~ d near the phase transition. For d = 3 this means M ~ t v / 2 . Note 
that this is exactly the scaling behavior of the true photon mass near the neutral 3D XY 
fixed point! This means that M can be alternatively regarded as a photon mass describing 
the crossover regime observed experimentally. This is the case if we interpret the GL model 
given in ( |II.1| ) as a disorder field theory that is, as a continuum dual GL model [fFHl2|. In 
this case, the charge e should be replaced by 2nM/q, where q is a the dual charge satisfying 
the Dirac condition qe = 2n. 

Let us investigate now how the exponent (3 should be defined. This presupposes a 
definition of order parameter whose scaling has exponent (3. This order parameter arises 
naturally from the definition of the superfluid density p s . This is defined for t < by 

Ps =< |0| 2 >= Z^ 1 < \<p \ 2 > . (11.26) 

We define the order parameter by $ by 



$ = J< |0 O | 2 >• (IL27) 



The expectation value of a gauge invariant operator is gauge independent p]20||. Thus, 
p s should be independent of a. However, is gauge dependent (see section III) and we 
have that $ should in fact depends on a to cancels the gauge dependence coming from Z^ 1 
in ( |11.26|) , otherwise p s would be gauge dependent. Note that the gauge invariance of 1 - 12 



implies that $ is independent of a but not necessarilly independent of a. However, its scaling 
behavior near the phase transition will be shown (section III) to be gauge independent or, 
more precisely, to be evaluated in the Landau gauge. The same will shown to be true for 
Zfr. For the moment, let us check that the above definition of order parameter works. The 
exponent (3 is defined through its behavior near the critical point, $ ~ \t\P . Also, we have 
\t\ vr) . Putting all of this together we obtain p s ~ \t\ 2 P~ UTJ , which is the Josephson 



m ' 



relation as obtained originally in Josephson's paper fl8 |. Using the hyperscaling relation 
dv = 2 — a together with the combination of the scaling relations a + 2/3 + 7 = 2 and 
7 = i/(2 — 77) to eliminate 7, we obtain the Josephson relation in the form p s ~ \t\ u ( d ~ 2 \ 



This last form of the Josephson relation has been also obtained directly by us [12], without 
using the scaling relations a + 2(3 + 7 = 2 and 7 = u(2 — 77). Therefore, 

(3= V -(d-2 + rj). (11.28) 

Thus, a gauge dependence in r\ would imply a gauge dependence in (3 if v is gauge indepen- 
dent. Note that the gauge independence of p s does not ensure the gauge independence of v. 
The point is that the gauge independence of \t\'^ d ~ 2:) can be a result of compensating gauge 
dependences arising from \t\ (equivalently, m§) and v. 



III. THE GAUGE DEPENDENCE 

In the preceding section we discussed the critical behavior of the GL model while trying 
to get some insight on the effect of the gauge dependence. Clearly a more careful analysis is 
needed. In this section gauge dependence always means a dependence on a, the renormalized 
gauge fixing parameter. Therefore, the other renormalized parameters of the model are 



7 



trivially gauge independent. This statement implies the gauge independence of all beta 
functions. Note that this does not means necessarilly that the bare parameters are gauge 
independent. 

In order to obtain the gauge dependence of we will employ the WT identities. From 
the WT identity Eq. ( |11.4| ) we obtain the following identity: 



M 



-d 2 z ) d*w!?\z;y,x) = ie \S(y - z)W (2) {z,x) - 6{z - x)W (2) {y, z) 



W {2 \p + k)-W {2 \p) 



By using twice Eq. (|III. 1|) we obtain 



W {2) 



d d k 



(k 2 + aM 2 y 



(III.l) 



(III.2) 



where W$ A ^(p) is the Fourier transform of 



W, 



(2) 



d d z 



< (d.ArfizWx^iy) > 



< (d,A,) 2 (z) >< <f>(x)4>\y) > 



(III.3) 



Let us denote the bare counterpart of W^ A ^ 2 by W^ A o^2 - 

,(2) 



We have that 



<9W (2) 



da 



(III.4) 



(9) i 

where W (x,y) =< (f>o{x)(j) (y) > is the bare 2-point connected correlation function. Eq. 
( pi.2|) is valid also if we replace the renormalized correlation functions by the bare ones and 
the renormalized couplings by their bare counterparts. Using then a bare version of (|III.2| ) 
and Eq. ( [111.4| ), we obtain 



dW[_ 
da 



(2) 



(P) 



d 3 k wi 2 \ P + k) - W^(p) 



(2), 



(27T) 



Eq. (pil.5|) can be rewritten as 



^W {2 \p) + — (p) 



da 

out of which we obtain 



da 



(k 2 + a Q M 2 ) 2 



d 3 k W ( - 2 \p + k)-W^ 2 \p) 
(2^p (P + a M 2 ) 2 ' 



(III.5) 



(III.6) 



d\nZ 4 
da 



-e n 



d d k 



1 



{2ir) d (k 2 + a M 2 ) 2 
d 



- 1 c d , 



(III.7) 



where n/Cd = (4ir) d / 2 r(d/2) sin(7r<i/2). Note that we have a pole for d — 4 in the second 
line of Eq. (|III.7| ). This is a consequence of the logarithmic divergence for d = 4. In the 
e-expansion the singular part of the different correlation functions is isolated as poles in 1/e 



S 



with e = 4 — d and the renormalization constants are written as power series in 1/e. This way 
of doing the things leads to the determination of the critical exponents as power series in e 
|J. The physical case of interest in critical phenomena of superffuid and magnetic systems 
corresponds to e = 1. 

By integrating the first line of ( |111.7| ) we obtain 

ln^(a ) = ln^(a = 0) - e 2 a / ^ y fc2(jfc2 -L^ . (III.8) 



Since egflo = e 2 a and clqMq = aM 2 , we can rewrite Eq. (|III.8|) as 

\nZ,{a) = lnZ,(a = 0) - e 2 a j + a M 2 ) 2 ' (IIL9) 

Let us assume that Z^a = 0) has been evaluated as a power series in 1/e. After regularizing 
dimensionally the integral in ( |111.9| ), we obtain 

% (a)=^(a = 0)-^, (111.10) 

which gives the gauge dependence of r\$ in the framework of the e-expansion. Let us assume 
that an infrared stable fixed point has been obtained, for instance, by ressummation methods 
||. As m — > 0, / — > /* 7^ (if e = 1), but a scales as m _1 near the fixed point and any 
non-zero a runs away as m — > 0. Thus, the only safe way towards the charged fixed point is 
over the line a = 0, that is, the Landau gauge. Note that for the case of interest in particle 
physics, d = 4, we obtain the same equation as ( |111.10| ). However, for d = 4 any gauge choice 
is possible since /* = in this case. 

In the fixed dimension approach things work differently. For d = 3 the integral in Eq. 
( |III.9|) is convergent and we can interchange the differentiation with respect to m with the 
integral sign. Since ae 2 and aM 2 are both RG invariants, we obtain 

^(a) = ^(a = 0), (Hill) 

and we obtain again that the physical gauge corresponds to a = 0. At this point some 
remarks are in order. First, from Eq. ( |111.11| ) we obtain drj^/da = while the same is not 
true for the rj^a) given in Eq. ( |III.10|) . Second, Eq. ( |IH.11| ) can be easily checked at 1-loop 
order. The renormalization constant is given as a function of /, a and v and if we take 
care of differentiating a when obtaining rj^, the result (|111.11| ) follows and coincides with Eq. 
( |II.21| ). Concerning the 1-loop example, it is instructive to ask ourselves what happens in 
other fixed dimension approaches. For instance, we could perform a critical point (m = 0) 
calculation where the renormalization conditions are defined at non-zero external momenta, 
taking the symmetrical point for functions which depend on more than one momentum 
variable [0,n. In this case the photon mass M is unecessary since the non-zero external 



momenta take care of infrared divergences [19]. The 1-loop expression for Za, in an arbitrary 



gauge is in this case rather simple and has been calculated by Schakel [21]. It turns out in 
this case that is independent of a if d = 3. 

The gauge dependence of Z m can be obtained in an analogous way. From the renormal- 
ization condition ( |IL8|) and m 2 = Z^Z^m^ we obtain that Wq (0) = Z m /m^. Using again 



9 



the bare version of (|III.2|) , we obtain exactly the same equation as Eq. ( |III.7|) but with Z$ 
replaced by Z m . This means that the gauge dependence of Z m is the same as for Z^. If we 
use the e-expansion we have that r\ m — r/^ is gauge independent since the gauge dependence 
of r] m will cancel exactly the gauge dependence of 7/0. In fixed dimension d = 3, on the other 
hand, r/ m (a) = r/ m (a = 0). It follows that the critical exponent v is gauge independent. Since 
i] is gauge independent, it follows that (3 is gauge independent and the order parameter $ 
has a true physical meaning. 



IV. CONCLUSIONS 



In this paper we have shown that the critical exponent (3 of the superconducting order 
parameter is gauge independent or, stating more correctly, that it must be evaluated in the 
Landau gauge a = 0. We may wonder if it is not a wasting of time to prove a result about 
something unaccessible experimentally. The point is that it is not sure that $ cannot be 
measured and we hope that the discussion in this paper could stimulate some experimental 
effort in this sense. Moreover, we have shown that transversality (the Landau gauge) is 
an intrinsic physical feature of the d = 3 GL model, a property not shared by the d = 4 
model. In d = 4 (the case of interest in particle physics) the Landau gauge is used due 
to its computational simplicity [E2]. In contrast, the Landau gauge is the only physically 
meaningful gauge in d = 3 [^3| . 
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FIG. 1. Flow diagram in the /a-plane. 
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